Operator splitting and approximate factorization for taxis-diffusion-reaction models
نویسندگان
چکیده
In this paper we consider the numerical solution of 2D systems of certain types of taxis-diiusion-reaction equations from mathematical biology. By spatial discretization these PDE systems are approximated by systems of positive, nonlinear ODEs (Method of Lines). The aim of this paper is to examine the numerical integration of these ODE systems for low to moderate accuracy by means of splitting techniques. An important consideration is maintenance of positivity. We apply operator splitting and approximate matrix factorization using low order explicit Runge-Kutta methods and linearly implicit Runge-Kutta-Rosenbrock methods. As a reference method the general purpose solver VODPK is applied.
منابع مشابه
Testing Weighted Splitting Schemes on a One-Column Transport-Chemistry Model
In many transport-chemistry models, a huge system of ODE's of the advection-diffusion-reaction type has to be integrated in time. Typically, this is done with the help of operator splitting. Operator splitting is attractive for complex large-scale transport-chemistry models because it allows to handle different processes separately in different parts of the computer program. Rosenbrock schemes ...
متن کاملOn the Efficient Treatment of Vertical Mixing and Chemistry in Air Pollution Modelling
Reaction-Advection-Diffusion problems occurring in Air Pollution Models have very stiff reaction terms and moderately stiff vertical mixing (diffusion and cloud transport) terms which both should be integrated implicitly in time. Standard implicit time stepping is by far too expensive for this type of problems whereas widely used splitting techniques may lead to unacceptably large errors for fa...
متن کاملNumerical Time Integration for Air Pollution Models
Due to the large number of chemical species and the three space dimensions, o -the-shelf sti ODE integrators are not feasible for the numerical time integration of sti systems of advection-diffusion-reaction equations @c @t +r (uc) = r (Kr c) +R (c) ; c = c(x; t); c 2 IRm; x 2 IR3 from the eld of air pollution modelling. This has led to the use of special time integration techniques. This paper...
متن کاملRobust numerical methods for taxis-diffusion-reaction systems: Applications to biomedical problems
In this paper we consider the numerical solution of taxis–diffusion–reaction models. Such nonlinear partial differential equation models appear often in mathematical biology and we consider examples from tumour growth and invasion, and the aggregation of amoebae. These examples are characterised by a taxis term which is large in magnitude compared with the diffusion term and this can cause nume...
متن کاملLinearly-implicit, approximate factorization, exponential methods for multi-dimensional reaction-diffusion equations
Linearly-implicit, approximate factorization methods are used to reduce the solution of multi-dimensional reaction–diffusion equations to a sequence of two-point boundaryvalue problems in ordinary differential equations which are solved analytically in a piecewise manner. It is shown that, depending on the allocation of the transient, diffusion and reaction terms to the second-order ordinary di...
متن کامل